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The phase behavior of ionic fluids on simple cubic and tetragonal (anisotropic) 
lattices has been studied by grand canonical Monte Carlo simulations. Systems with 
both the true lattice Coulombic potential and continuous-space 1/r electrostatic in- 
teractions have been investigated. At all degrees of anisotropy, only coexistence 
between a disordered low-density phase and an ordered high-density phase with 
the structure similar to ionic crystal was found, in contrast to recent theoretical 
predictions. Tricritical parameters were determined to be monotonously increasing 
functions of anisotropy parameters which is consistent with theoretical calculations 
based on the Debye-Hiickel approach. At large anisotropics a two-dimensional-like 
behavior is observed, from which we estimated the dimensionless tricritical temper- 
ature and density for the two-dimensional square lattice electrolyte to be T^* ^ = 0.14 
and pl^i = 0.70. 
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I. INTRODUCTION 

After controversial experimental results on the nature of ionic criticality the 
studies of critical phenomena in ionic fluids received a fresh impetus, and over the last 
decade substantial progress has been made However, a complete thermodynamic 

description of electrolyte systems is still unavailable. 

The simplest and most frequently employed model system for both theoretical and Monte 
Carlo investigations of Coulombic systems is the restricted primitive model (RPM), which is 
a gas of equisized hard-sphere ions carrying positive and negative charges of equal magnitude. 
From the theoretical point of view, a reasonable description of the RPM in the critical region 
has been obtained ^t a mean- field level using integral equations approach p,l2l, 



Hiickel theory 



and Debye- 
911 agree 



While theoretical predictions of the critical parameters 
reasonably well with those obtained in computer simulations |lo[ recent Monte Carlo 
investigations of ionic systems indicate a non-classical Ising-like asymptotic critical behavior 
jl^ . Moreover, predictions of different mean-field theories for charge- and size-asymmetric 
primitive systems differ significantly from each other Q, Q, Q, In the absence 
of real experimental data, numerical simulations provide an important criterion for testing 



different theoretical methods |l8, 



y 



2l|. 



In recent years, a new theoretical approach to study the thermodynamics and criticality 
of Coulomb systems based on the investigation of corresponding lattice models, has been 
introduced [l^, 22, 2^, 24, Being obviously less realistic than their continuous-space 

counterparts, lattice models are much easier to handle analytically, and the information they 
provide is crucial for understanding the physics of criticality. Moreover, finely discretized 
lattices with lattice spacing several times smaller than the ion size present substantial com- 
putational advantages over continuous space simulations, while producing almost identical 
critical parameters and phase diagrams 1^. By contrast, when the ionic diameter is equal 
to the spacing of the underlying lattice, the phase diagram is drastically different. The sym- 
metry of simple cubic or body-centered cubic lattices allows for unequal charge distribution 
between two identical sublattices, yielding an ordered structure similar to an ionic crystal. 
The competition between electrostatic interactions of the two charged sublattices and the 
entropy of charge separation leads then to order-disorder phase transitions and a tricritical 



point 



mm 
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A systematic study of the lattice restricted primitive model (LRPM) based on the Debye- 
Hiickel approach has shown that, in accordance with Monte Carlo simulations, a Coulombic- 
governed gas-liquid phase transitions seen in continuum electrolytes is totally suppressed on 
simple cubic and body-centered cubic lattices due to the formation of the thermodynamically 

n n 

more stable charge-ordered phase |2j]. Ciach and Stell |23] argued that, unlike non-ionic 
fluids, in charged systems the most important are short-range fluctuations. Consequently, 
they studied a LRPM supplemented with additional short-range interactions between ions 
and concluded that at some strength of the short-range attraction both the normal critical 
and tricritical point can become thermodynamically stable 2^ • This leads to complex phase 
diagrams with two distinct phase transitions, one between gas and disordered liquid phase, 
and another one between disordered liquid and a charge ordered phase. Recent Monte 
Carlo studies have confirmed qualitatively most of these theoretical predictions but further 
investigations are required j^]. 

At low dimensions Coulombic interactions are stronger than in three-dimensional systems, 
and the possibility for gas-liquid phase transitions increases. This idea has been used recently 
in the investigation of anisotropic lattice models of electrolytes using Debye-Hiickel method 
j^l • In this study the anisotropy mimicked the lowering of the spatial dimensionality. It has 
been found, that for strongly anisotropic lattices gas-liquid phase coexistence is restored. 
However, this theory neglects an ion clustering phenomena, which are important in charged 
particles systems 0, 0]. Thus these theoretical predictions should be tested by computer 
simulations. 

Current Monte Carlo studies of lattice RPM in ci = 3 dimensions employ the contin- 
uous space 1/r potential HHQ- 

However, because of the discrete lattice symmetry, 
the correct lattice Coulomb potential differs from 1/r at short distances and approaches it 
asymptotically at large distances. This raises the question of how using the correct elec- 
trostatic interactions will affect the thermodynamics and critical behavior of ionic systems. 
Indeed, the lattice correction to 1/r potential decays rapidly with distance and thus it can be 
viewed in a sense as an additional short-range interaction. In order to understand correctly 
the underlying mechanisms of phase transitions in Coulombic systems, lattice models should 
be investigated with the corresponding lattice potentials. 

In this paper we study the RPM on tetragonal lattices using grand canonical Monte 
Carlo simulations and histogram reweighting technique [s^. The correct lattice Coulomb 
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potential is used and an Ewald-like summation is utilized to account for the long-range 
nature of electrostatic interactions. The paper is organized in the following way. In Sec. 
II the details of the model and simulation method are given. The results and discussions 
are presented in Sec. Ill, and summary and conclusions are drawn in Sec. IV. Analytical 
expressions and the details of lattice potential calculations are given in Appendix. 



II. ANISOTROPIC LATTICE MODEL AND SIMULATIONS METHOD 



The system we have studied consists of 2N ions, half of them carrying charge q and 
half — g, positioned on a simple tetragonal lattice with lattice parameters ratio (degree of 
anisotropy) a = hja. The pairwise interaction potential of two ions separated by distance 
Tij has the form (for details see Appendix A) 



An'^Db J-TT 2 + l/a^ — (cos ki + cos k2 + l/a"^ cos 

where D is the dielectric constant of the structureless solvent. Reduced quantities are 
obtained by scaling with the energy of the strongest ion-ion interaction in two neighboring 
xy planes in the continuum limit, Eq = q'^/Db, and the unit cell volume, vq = o?h, namely, 

where V is the volume of the system. 

Since the Coulomb potential is long ranged, one needs to account for interactions with 
particles in all images of the periodic box used in the simulations. In the case of continuous- 
space 1/r interactions the standard method to achieve this is the Ewald summation (see 



e.g. Ref. 



281 



29|). However, no similar method is available for the lattice Coulomb potential. 



Therefore, we propose the following approximate scheme, which is based on the asymptotic 
properties of the anisotropic lattice potential ((T)). 

It can be shown that the potential ((H) for finite a at large distances behaves as 

^.,(r=(x,,,.))~-^=^=^_, (3) 

For this potential, an Ewald-like sum can be constructed along the lines of the derivation in 



Ref. 



29l | , which yields the "continuum" part of the energy of interactions of periodic system 
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Table I: Dimensionless lattice Coulomb potential Uiatuce © and its asymptotic limit Ucont given 
by Eq.® 



cry st allogr aphic 


a- 


=1 




=4 




=10 


indexes 


Ulattice 


Ucont 


Ulattice 


Ucont 


Ulattice 


Ucont 


[10 0] 


1.0815164 


1.0000000 


3.192933 


4.000000 


4.949445 


10.00000 


[0 1] 


1.0815164 


1.0000000 


1.028176 


1.000000 


1.006922 


1.000000 


[10 0] 


0.1002578 


0.1000000 


0.410816 


0.400000 


1.066964 


1.000000 


[0 10] 


0.1002578 


0.1000000 


1.000167 


0.100000 


0.100015 


0.100000 


[10 10 10] 


0.0577029 


0.0577350 


0.094250 


0.094281 


0.099021 


0.099015 


[100 0] 


0.0100003 


0.0100000 


0.040009 


0.040000 


0.100129 


0.100000 


[0 100] 


0.0100003 


0.0100000 


0.010003 


0.010000 


0.010041 


0.010000 


[100 100 100] 


0.0057736 


0.0057733 


0.009431 


0.009428 


0.009924 


0.009901 



of 2N particles 

^ ^ 1 ^ f^a^exp{-n^\n'\yri^ + 2mn'v[A 1 „ erfc(r/lr,, + n'] 



all space I ijl «J=1 \ n' II n' I *J I / 

where r' = {x/a,y/a,z), n' = {anx,any,nz), < nx,ny,nz < 27r, and erfc(x) is the 
comphmentary error function. Here conducting boundary conditions have been utihzed, 
which are less sensitive to finite-size effects |19]. As usual, the real-space damping parameter 
7] must be chosen in a such way that all sums in (0} converge fast. In our simulations we 
used 1] = 5, and the Fourier-space sums were restricted to 518 wave vectors. 

The difference between the correct lattice potential and its asymptotic limit decays 
rapidly with distance: see Table 1. Therefore, to obtain the total energy of electrostatic 
interactions, the lattice correction is added for ions in elliptic shells up to certain distance, 

Etotai = E,+ y: ff^^^W - t;^) • (5) 

\r'/L\<nraax ^ I «i ' / 

We performed grand canonical Monte Carlo simulations on cubic boxes of length L 
stretched in z direction by a factor of a, under periodic boundary conditions. Distance- 
biased algorithm was employed for insertions and removals of pairs of unlike ions at each 
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time step in order to facilitate acceptance, following Ref.j31|. To analyze the simulation 
data and obtain the coexistence curves, multi-histogram reweighting techniques js^ have 
been used. 

III. RESULTS AND DISCUSSION 

To reduce possible finite-size effects, we used a system size L* = L^^y/a = L^/b = 12 
for the isotropic lattice {a = 1) and a = 1.2, and L* = 16 for stronger anisotropics. Short 
initial runs were used to get a rough estimate of the location of the tricritical point. After 
that, typical runs involved 1 — 3 x 10'' Monte Carlo steps for equilibration and 2 — 7 x 
10^ for production. We did not try to locate the second order transitions Neel line and 
restricted ourselves with obtaining approximate location of the tricritical point by simple 
linear extrapolation of the coexistence lines, as expected for d = 3 tricriticality. 

In order to separate the effect of using the correct lattice potential in simulations, we 
studied first the phase coexistence in LRPM on a simple cubic lattice: see Fig.l. On a 
simple cubic lattice. Coulomb gas with 1/r potential is known to phase separate, with two 
coexisting phases being a low-density disordered phase and an ordered high-density ionic 
crystal-like phase Q,|2^. The tricritical par ameters estimates for the system size L* = 12, 
obtained by Panagiotopoulos and Kumar [lOj , are T*^^ = 

0,15±0,01.p;„ = 0.48 ±0,02 0. 



For the same system size, taking into account the lattice correction to the potential yields 
T^* j = 0.22, p^^j = 0.48 with the same accuracy. Thus the tricritical density remains the same 
within the error limits, but the temperature increases by the factor of 1.5. Qualitatively, the 
increase in the tricritical temperature should be expected since lattice Coulomb interactions 
are stronger than 1/r at short distances, which leads to a higher stability of dense phases. 
Nevertheless, the tricritical point is still significantly lower ( by 5 0%) than the predictions 
of the Debye-Hiickel theory of lattice electrolytes Ttri ~ 0.365 [2^, but is about 20% higher 
than the results Tt^i ~ 0.202 from the hierarchical reference theory (HRT) of LRPM by 
Brognara et al.j2^, who also used the correct form of potential in their calculations. Note, 
that both of these theories agree on the tricritical density to be p^^.^ ~ 0.38, which is about 
20% lower than our simulation results. However, while taking into account Bjerrum pairing 
of ions into neutral dipoles |3^ in Debye-Hiickel approach would supposedly produce better 
results by shifting the phase coexistence to lower temperatures and higher densities, the 
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predictions of the HRT seem to be final. 

The phase coexistence for different degrees of anisotropy a is shown in Fig.l. At any 
a, the system shows only order- disorder phase separations and a tricritical point. Thus, 
the effect of stretching the lattice amounts merely to altering the tricritical point location, 
without changing qualitatively the topology of the phase diagram. This goes in contrast with 
the recent theory based on Debye-Hiickel calculations 26] which predicts that at the lattice 
parameters ratio h/a> 2.6 a distinct gas-liquid phase coexistence may reappear. However, 
as was mentioned above, the theoretical approach of Ref.j^l disregards the Bjerrum pairing 
or, more generally, clustering of ions, which is crucial for the thermodynamics of the gas- 
liquid transition in electrolytes. Typically, taking into account pairing somewhat diminishes 
the critical temperature but significantly increases 

taking into consideration the ion clustering is unlikely to change the topology of the phase 
diagram of LRPM on the isotropic lattice. However, since the strength of the electrostatic 
interactions increases with anisotropy, this may lead to more profound tendency of ions to 
bind into neutral dipoles/clusters (mostly, in xy-planes), producing eventually a Kosterlitz- 
Thouless (KT) transitions line at infinite stretching. Stronger pairing then could shift 
the gas-liquid coexistence curve much further than in the isotropic case, again putting it 
completely inside the order-disorder phase envelope. It should be noted that simulations 
of the two-dimensional lattice Coulomb gas 3^ also reveal no normal first order gas-liquid 
transition, but only an order-disorder phase separations [s^. 

Tricritical temperature and density as functions of the lattice parameters ratio a are 
presented in Fig. 2. As anisotropy increases, both tricritical temperature and density in- 
crease. Physical account of this behavior follows from the properties of the lattice Coulomb 
potential. Indeed, in z-direction lattice potential behaves almost as I/2; (see Table 1). 
However, in crystallographic directions different from [001], the anisotropic lattice potential 
is much stronger than 1/r and decays slower. While the interactions in 2;-direction almost do 
not change with a, stronger anisotropy brings stronger interactions inside the xy-layers. This 
becomes especially important at large densities, where the average distance between ions is 
of the order of the lattice spacing. Therefore, the average energy per particle increases with 
anisotropy. As a consequence, higher temperatures are required for fluctuations to achieve a 
critical level, and the tricritical temperature grows. At the same time, more intensive inter- 
actions between ions lead to more active clustering. The free charges density, and hence the 
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electrostatic part of the free energy diminishes significantly. However, the entropic part of 
the energy remains the same. Therefore, electrostatic interactions of two charged sublattices 
become comparable with the entropy of charge separation at higher overall densities, and 
the coexistence shifts to denser phases. 

The trend for tricritical temperature and density to increase with anisotropy is correctly 
reproduced by the Debye-Hiickel treatment of anisotropic lattices 3|: see Fig. 2. However, 
the quantitative agreement between theoretical predictions and simulations is not very good. 
While the theory predicts only a few percent increase in density at large lattice stretching, 
our simulations showed more significant growth. For the tricritical temperature, on the 
contrary, simulations yield a somewhat smaller slope. 

At large stretching, increasing interactions inside xy-planes result in lesser relative impor- 
tance of interactions between ions with different z coordinates. At some point, interactions 
with ions in the same layer start dominating, and the 3D system becomes quasi-2D. Since 
different layers become uncoupled, due to its probabilistic nature phase separation in each 
of them takes place at different values of the chemical potential, and one may have config- 
urations with gas phase in some of the layers and high-density ordered phase in the others, 
which prevents reliable sampling of the system. At a=10, it turned to be impossible to 
obtain a reliable phase separation in a periodic system of more than one layer in ^-direction: 
at high densities the system would not equilibrate after more than 10^ Monte Carlo steps. 

To this end, for a = 10 we carried out simulations of the system consisting of only 
one plane, that is, square box with = L* = 32 periodic in x and y directions, with 
the potential (^. Using only one layer in simulations yielded a highly asymmetric phase 
coexistence very similar to that obtained by Teitel for Coulomb gas on two-dimensional 
square lattice 35[, with the tricritical point located at T^*- = 0.28, pj'^j = 0.70: see Fig.l. 
But in two dimensions, the coefficient in the Laplace equation ()A1|) is not Air but 2tt. Since 
all our simulations were done with three-dimensional potential, the temperature at large 
anisotropics must be divided by 2 to obtain the correct 2D limit, while the density should 
remain intact. This yields = 0.14, = O-'^O which exactly coincides with the values 
obtained by Teitel. 

It should be also noted that parts of the coexistence curve for a = 10 at intermediate 
densities lay below the corresponding curve for a = 4. However, for the latter case in contrast 
to a = 10, the interactions between layers are still very important and add significantly to 
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the total free energy of the system: simulations of periodic system of only one layer produce 
much lower coexistence curve. Therefore, although at equal densities the energy of ions 
interactions inside a single layer for a = 4 is less than for a = 10, the total energy density 
may be higher because of the interactions between layers, and hence the temperature of 
phase separations would also be higher. However, this is not essential at very low or high 
densities where the interactions inside layers dominate anyway. Taking into account the 
influence of layers on each other at a = 10 probably would eliminate this intersection of the 
coexistence curves. Nevertheless, this hardly would change our estimate for the tricritical 
point location in the limit of infinite stretching, which, as it has been mentioned, coincides 
with the results by Teitel for two-dimensional square lattice. 

To gain further insight on the behavior of LRPM on anisotropic lattices we also per- 
formed simulations with the continuous space Coulomb potential. The corresponding phase 
diagrams for box size L* — 16 and a > 2 are presented in Fig. 3. At any fixed lattice stretch- 
ing, the phase coexistence is very similar to the case of the true lattice potential. No normal 
gas-liquid phase transition is found, and the phase coexistence is between a disordered low- 
density phase and ordered high-density phase. The tricritical temperature shows the same 
trend as for the correct lattice potential, increasing with the inter-layer distance. However, 
since for small a the interactions are weaker now, the tricritical temperatures are lower. 

The scaling of the temperature T with Eq — /Dh is maintained in this paper to allow 
comparisons with earlier theoretical results. For the correct lattice (or continuum) two- 
dimensional potential, this is allowed as 6 — > oo because the potential becomes length scale- 
independent at the two-dimensional limit. However, for the 1/r potential, using a scaling 
with the layer-layer energy results in a divergence of the reduced critical temperature. While 
for the correct lattice potential the "natural" scale of energy is set by layer-layer interactions 
whose asymptotic behavior does not change with a (see Table 1), for the 1/r potential this 
scale is given by the energy of the ion- ion interactions at contact inside a layer Eq — / Da. 
Note, however, that for actual simulations this is not important, because all what matters 
is the Boltzmann factor, which remains the same as long as the energy and the temperature 
are measured in the same units. 

Calculations with the temperature scaled by this energy T" — kTDa/q^ (where D is the 
three-dimensional dielectric constant, appropriate for the 1/r potential) are shown in Fig. 
4. With this "natural" scaling, the tricritical point behavior as a function of anisotropy 
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is drastically different: the tricritical temperature decreases with anisotropy, converging 
to T" = 0.11. The reason for this is that now increasing anisotropy is not accompanied by 
stronger interactions between ions. On the contrary, while interactions inside the xy— planes 
remain the same, due to a weaker coupling of the layers effective energy per particle dimin- 
ishes, and so does the tricritical temperature. This decay is rather fast and the temperature 
saturates already at the lattice parameters ratio o; ~ 2. Detailed studies of ions configura- 
tions showed that at a > 2 different xy— planes become totally independent of each other, 
and the system behaves as a quasi-2£) liquid. Similarly to the true lattice potential case, to 
obtain phase coexistence for the limit of infinite stretching (a > 2), simulations of periodic 
system of only one plane were carried out with box size L = 32 and L = 48. 

The tricritical density gradually increases from p^,.,j ~ 0.48 for isotropic lattice to pj'^j ~ 
0.55 at strong anisotropy, following the trend shown by the system with correct lattice 
interactions: see Fig. 3. This again reflects the relative importance of clustering when one 
approaches the two-dimensional limit of infinite stretching. However, since the potential is 
now weaker, the limiting density is lower than in true lattice 2D systems. 



IV. CONCLUSIONS 



In this study, grand canonical Monte Carlo simulations and histogram reweighting tech- 
niques have been used to investigate phase behavior of the restricted primitive model with 
correct lattice Coulomb potential on simple cubic and tetragonal lattices. Our results show 
that, at all degrees of anisotropy, only order-disorder phase separations and a tricritical 
point exist. For isotropic cubic lattice, the effect of using the correct lattice potential shows 
only in the value of the tricritical temperature. As anisotropy increases, both the tricritical 
temperature and density increase. At large anisotropy the system undergoes a qualitative 
shift from a three-dimensional to two-dimensional behavior. This shift is not accompanied 
by any change in the phase diagram topology, and only the tricritical parameters alter. To 
study in detail the transition from 3D to 2D behavior, further investigations are needed. 
In particular, its location may depend on the system size and boundary conditions used 
in Ewald summation. For continous 1/r potential, the tricritical temperature displays an 
opposite trend and decreases with the anisotropy. Although exact values of the tricritical 
temperature and density can be influenced by finite size effects, an issue totally omitted in 
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this study, their quahtative behavior seems physically plausible. It would be also useful to 
develop another procedure for Ewald-like summation for lattice potential, which will not not 
rely on its asymptotic behavior. 

The only existing theoretical treatment of Coulomb gas on anisotropic lattices based 
on the Debye-Hiickel approach captures ^ correctly the trend of the tricritical temperature 
and density behavior with anisotropy [26]. However, the absence of the anticipated normal 
gas-liquid transition suggests that with increasing anisotropy, stability of both ordered and 
disordered dense phases increases in a similar way, and the formation of an ordered struc- 
ture remains thermodynamically more favorable than possible gas-liquid phase coexistence. 
Taking explicitly into account the ion clustering will provide a better theoretical description 
of thermodynamics of ionic systems on anisotropic lattices. 
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Appendix A: LATTICE COULOMB POTENTIAL: EXACT REPRESENTATION 

AND FAST CALCULATION 

Due to the discrete symmetry, the correct lattice Coulomb potential differs from continous 
1/r. The analytical expression for this potential follows from the lattice version of the Poisson 
equation 23^ . 

A^{r) = -^5(r), (Al) 

where the exact form of the lattice Laplacian depends on the geometry of the lattice, D is 
the dielectric constant of the media, and Vq is the unit cell volume. For simple tetragonal 
lattice with lattice spacings a in x and y directions and b in z direction one has 

Aif = A^^ + Ay^ + A,^, (A2) 
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with 

Aiip{r) = \la\ [(^(r - 0^6^) - 1ip{v) + v?(r + a^ei)] , (A3) 

where i 6 and are the unit vectors along the corresponding 

lattice directions. After introducing the anisotropy parameter a = 6/a, the solution of the 
Poisson equation has the following form 

^ 1 r e^'^'-ci^k 

^^^^ ~ iTi^Db J -^2 + 1/ -{cosh + cos k2 + 1/0"^ cos ksY 

However, the potential in the form of a triple integral of a periodic function is not con- 
venient for numerical calculations. More practical expression for this integral has been 
obtained by Maradudin et al. 3^ 

V?(r = (/, m, n)) = exp[-(2 + l/a^)t]Ii{t)I^{t)In{t/a^)dt, (A5) 

where Ii denotes a modified Bessel function of the first kind. 

To enhance calculations of the lattice potential for the whole lattice, an error-free propa- 
gating algorithm, similar to that proposed by Friedberg and Martin for isotropic lattice, 
has been used. First, we compute lattice potentials in planes {xOz) and {xlz) using Eq. ()A5|l . 
Now, for a given lattice cite {k,l,m), the anisotropic lattice potential satisfies the gradient 
equations 

ip{l, m, n) — ipil — 1, m, n) ip{l, m + l,n) — ip{l, m — 1, n) 

/ m 

2(p{l,m,n + l) - (p{l,m,n-l) 

= a (Ad) 

n 

which follow from Eq. ()A4jl . At fixed n, potential at any point {l,m,n) can be obtained 
recursively from two lines, m = and m = 1. For / = it would be indeterminate, but 
this is not important since due to the symmetry (f{0,m,n) = (f{m,0,n). As Friedberg 
and Martin showed [s^, this recursion does not lead to increase of the initial possible 
rounding error at m = and m = 1 lines as long as |/| < |m|. For |Z| > |m| the symmetry 
ip{l, m, n) = ip{m, /, n) can be employed. 

Note that similar considerations can be done for initial (xO^;) and {xlz) planes in attempt 
to generate the potentials from two lines, using the second equation in ()A6|) . In these planes 
the potential is not symmetric, and therefore two sectors, |/| < a^\n\ and |/| < a^\n\ need to 
be calculated separately, starting from the corresponding axes and then " stitched" together 
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at the \l\ — a^\n\ line. However, since all potentials and interactions are computed only 
once at the beginning of simulations and then fast look-up algorithm is used, this gives us 
a relatively small advantage in the overall computational time. 
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Figure 1: Phase diagrams in LRPM with correct lattice potential on tetragonal lattices for different 
degrees of anisotropy (open circles). Open circles from bottom to top a = 1, 1.2, 2, 4; solid triangles 
correspond to a=10. For reference, phase coexistence in LRPM on cubic lattice with 1/r potential 
[lol | is shown in solid circles (see also Fig. 3). 
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Figure 2: Reduced tricritical parameters as a function of the lattice spacing ratio, (a) tricritical 
temperature, (b) tricritical density. Open triangles correspond to the correct lattice potential, 
squares are for the continuous space potential 1/r. Dashed lines are merely guides to the eye. The 
analytical theory predictions are shown in solid lines |2^. 
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Figure 3: Phase diagrams of LRPM with continuous-space 1/r Coulomb potential. Prom bottom 
to top, a = 1, 1.2, 1.5, 2. Reduced temperature is defined by Eq.(2). 
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Figure 4: Reduced tricritical temperature T" = kTDa/q^ for the system with continuum 1/r 
Coulomb potential as a function of the anisotropy. 



